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1 .  INTRODUCTION  AND  BACKGROUND 

The  Spaced  Based  Infrared  Systems  (SBIRS)  Program  is  a  multi-satellite  constellation  designed  to 
conduct  global  and  regional  surveillance  using  passive  infrared  detection  techniques.  To  facilitate  the  smooth 
operation  of  instruments  on  board  satellites,  the  impacts  of  clouds  as  a  clutter  source  and  as  obscurants  must 
be  quantified.  One  of  the  primary  goals  for  the  cloud  related  efforts  is  to  assess  the  spatial  variability  of 
background  radiances  in  spectral  bands  between  2  and  5  pm  at  geometries  relevant  to  the  SBIRS  constellation. 
Clouds  are  important  contributors  to  the  overall  background  radiance  structure,  particularly  in  spectral  bands 
with  weighting  functions  in  the  troposphere  and  lower  stratosphere.  In  particular,  the  high-altitude  cirrus 
clouds  often  pose  a  serious  problem  in  the  satellite  surveillance  using  the  passive  infrared  detection  technique. 
In  this  project,  we  have  carried  out  a  program  designed  to  understand  the  effects  of  cirrus  clouds  between  2 
and  5  pm  on  the  upwelling  radiances  measured  by  infrared  sensors  on  board  satellites,  and  to  develop  a  cirrus 
cloud  remote  sensing  algorithm  based  on  their  physical  properties  within  this  spectral  range  for  application  to 
the  development  of  infrared  sensors. 

The  principal  objectives  of  the  present  project  are:  (a)  to  compute  the  ice  crystal  phase  function, 
single-scattering  albedo,  and  extinction  coefficient  for  the  spectral  range  2-6.4  pm  using  a  unified  theory  for 
light  scattering  by  ice  crystals  developed  by  our  research  group;  (b)  to  collect  and  analyze  available  ice  crystal 
size  distributions  to  support  the  radiance  measurements  from  satellite  and  aircraft;  (c)  to  develop  a  remote 
sensing  algorithm  and  computational  code  for  the  retrieval  of  cirrus  temperature  (height),  optical  depth,  and 
ice  crystal  size  based  on  the  radiative  properties  of  cirrus  clouds;  and  (d)  to  collaborate  with  the  SBIRS 
scientists  to  analyze  the  available  aircraft  and  nucrophysic  data.  To  meet  the  preceding  research  goals,  we  have 
acquired  the  radiance  data  collected  by  the  75 -channel  Airborne  Remote  Earth  Sensing  (ARES)  System  and 
the  microphysical  data  collected  by  the  cloud  probes  on  board  the  High  Altitude  Reconnaissance  Platform 
(HARP)  during  the  Hanscom  AFB  mission  flight  on  September  16,  1995. 

This  report  is  organized  as  follows.  In  Section  2,  we  discuss  the  range  of  ice  crystal  bounding 
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parameters,  including  the  refractive  indices  for  ice  in  the  ARES  spectral  interval  and  the  ice  crystal  size  and 
shape  for  midlatitude  cirrus  cloud  systems.  Analyses  of  the  available  ice  crystal  size  distribution  measurements 
were  made  to  obtain  representative  values  for  radiative  transfer  and  remote  sensing  applications. 
Computational  results  for  the  extinction  coefficient,  single-scattering  albedo,  asymmetry  factor,  and  phase 
function  for  the  ARES  75 -channel  spectral  wavelengths  and  for  ice  crystals  of  various  sizes  and  shapes  are  also 
presented.  Effects  of  incorporation  of  small  ice  crystals  on  scattering  and  absorption  calculations  are  further 
investigated.  Section  3  describes  the  development  of  a  remote  sensing  algorithm  for  the  retrieval  of  cirrus 
temperature  (height),  optical  depth,  and  ice  crystal  size  using  the  ARES  5. 1-5.3  and  3.7  pm  channels. 
Validation  of  the  retrieval  scheme  is  further  made  by  using  the  ground-based  8.6-mm  radar  and  aircraft  2D-C 
probe  ice  crystal  data.  Finally,  summary  is  presented  in  Section  4. 

2.  LIGHT  SCATTERING  BY  ICE  CRYSTALS 
2 . 1  Index  of  Refraction  for  Ice 

For  interpretation  of  the  observed  ARES  bidirectional  reflectances  (BRDFs)  and  development  of  a 
remote  sensing  algorithm  based  on  ARES  channels  to  retrieve  cirrus  cloud  temperature,  optical  depth,  and  ice 
crystal  sizes,  we  must  understand  the  spectral  variation  of  the  refractive  indices  for  ice,  the  ice  crystal  size 
distribution,  and  the  variability  of  ice  crystal  shape. 

Figure  1  shows  the  real  and  imaginary  refractive  indices  for  ice  and  water  from  2  to  6.5  pm  (Liou 
1992).  The  real  part  of  the  refractive  index  is  a  measure  of  the  phase  speed  of  the  electromagnetic  wave,  while 
the  imaginary  component  of  the  refractive  index  is  related  to  the  absorption  coefficient.  The  drastic  increase 
of  the  real  part  from  about  2.9  to  3.2  pm  for  both  ice  and  water  is  referred  to  as  the  anomalous  dispersion.  The 
imaginary  part  varies  by  about  three  orders  of  magnitude  from  2  to  3.5  pm.  For  A  >  3.5  pm,  the  variation  is 
smaller.  Both  water  and  ice  show  similar  patterns  of  the  spectral  variation.  To  resolve  the  aforementioned 
spectral  variation  of  the  refractive  indices  of  ice,  we  have  selected  the  optimal  number  of  wavelengths  for  light 
scattering  calculations  that  are  required  for  applications  to  the  ARES  channels.  More  wavelengths  are  chosen 
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for  the  spectral  region  where  the  refractive  index  varies  drastically.  Appropriate  interpolations  can  be 
performed  to  obtain  the  scattering  and  absorption  values  corresponding  exactly  to  the  ARES  channel 
wavelengths. 

2.2  Ice  Crystal  Size  Distributions 

The  interpretation  of  radiometric  measurements  and  the  development  of  remote  sensing  algorithms  for 
detecting  ubiquitous  cirrus  clouds  requires  comprehensive  knowledge  of  the  two  fundamental  ice  microphysics 
parameters:  ice  crystal  size  distribution  and  ice  crystal  shape.  Both  parameters  vary  significantly  in  space  and 
in  time  and  their  variabilities  are  subject  to  a  number  of  ice  microphysical  processes  spanning  from 
homogeneous  and  heterogeneous  nucleations  to  diffusional  growth,  collision  and  coalescence,  and  gravitational 
settling  processes.  Because  aircraft  and  sounding  measurements  are  restricted  in  their  spatial  and  temporal 
coverages,  only  limited  databases  are  available.  In  the  following,  we  present  three  representative  ice  crystal 
size  distributions  which  are  chosen  from  various  ice  crystal  size  distributions  collected  by  replicators  and  cloud 
probes  during  a  number  of  field  experiments  over  midlatitude  regions. 

In  conjunction  with  our  radiative  transfer  and  remote  sensing  studies  involving  cirrus  clouds,  we  have 
analyzed  11  ice  crystal  size  distributions  from  aircraft  observations  for  use  in  the  parameterization  of  their 
radiative  properties  (Fu  and  Liou  1993).  We  have  also  acquired  aircraft  optical  probe  data  and  balloon-borne 
replicator  sonde  data  collected  during  FIRE-II-IFO,  Kansas,  1991,  and  during  April  1994  over  the  ARM- 
CART  site.  From  these  ice  crystal  size  distributions,  we  have  selected  three  representative  cases:  the  balloon¬ 
sounding  data  collected  during  FIRE-II-IFO  on  December  5,  1991  and  two  ice  crystal  size  distributions 
corresponding  to  the  cold  and  warm  cirrus  clouds,  as  shown  in  Figure  2.  The  December  5  curve  is  the 
composite  ice  crystal  size  distribution  obtained  from  averaging  over  the  ice  crystal  size  distributions  at  selected 
vertical  levels  normalized  by  the  observed  number  density  at  each  level  (Ou  et  al.  1995).  On  that  date,  a  layer 
of  well-defined  cirrus  cloud  was  present  over  the  observation  site.  Near  the  cloud  top  at  13  km,  the  cirrus  cloud 
was  composed  of  small  ice  particles  with  sizes  between  5  and  100  pm.  Note  that  the  smallest  ice  crystals  that 
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Figure  2  Three  representative  ice  crystal  size  distributions  for  typical  midlatitude  cirrus  cloud  systems. 


the  replicator  could  resolve  are  on  the  order  of  5  pm.  In  the  middle  part  of  the  cloud  between  10  and  12  km, 
ice  crystals  were  columns,  bullet  rosettes,  plates,  and  aggregates.  The  size  distributions  become  broadened  and 
show  a  bimodal  shape  toward  the  lower  part  of  the  cloud.  Near  the  cloud  base  at  9.6  km,  the  size  distribution 
is  dominated  by  partially  sublimated  irregularly  shaped  ice  crystals.  The  composite  ice  crystal  size  distribution 
spans  from  about  5  to  800  pm,  which  is  the  typical  range  of  ice  crystal  size  distributions  observed  in  cirrus 
clouds.  Because  of  the  typical  size  range  and  the  presence  of  the  frequently  observed  ice  crystal  shapes,  the 
December  5  case  is  selected  as  the  representative  ice  crystal  size  distribution  for  midlatitudes. 

The  ice  crystal  distributions  associated  with  cold  cirrus  were  obtained  from  the  FSSP  measurements 
over  the  ARM-CART  site  at  about  0236  UTC,  April  19,  1994.  The  image  of  the  collocated  polarization  lidar 
return  signals  (K.  Sassen,  private  communication)  indicated  that  a  layer  of  cirrus  cloud  was  located  between 
12  and  13.5  km.  The  University  of  North  Dakota's  Citation  aircraft  was  flying  at  13  km,  where  the  measured 
temperature  was  -70°C.  The  cold  cirrus  size  distribution  extends  from  about  5  to  55  pm  and  fits  well  to  a 
Gamma  distribution.  The  ice  crystal  size  distribution  associated  with  warm  cirrus  clouds  is  taken  from 
Heymsfield  and  Platt  (1984)  which  was  derived  from  die  measurements  of  FSSP,  1D-C,  and  1D-P  probes  with 
temperatures  in  the  range  between  _25  and  -20°C.  The  size  ranges  from  about  20  to  1800  pm.  The  preceding 
three  ice  crystal  size  distributions  are  used  for  the  computations  of  the  single-scattering  properties  of  cirrus 
clouds  associated  with  the  ARES  channels. 

2.3  Scattering  and  Absorption  Properties  of  Ice  Crystals 

An  exact  solution  for  the  scattering  of  light  by  hexagonal  ice  crystals  covering  all  shapes  and  sizes  does 
not  exist  in  practical  terms.  Several  numerical  methods  that  have  been  developed  to  solve  the  problem  of  light 
scattering  by  nonspherical  particles  are  usually  applicable  to  size  parameters  smaller  than  approximately  20, 
as  discussed  by  Liou  and  Takano  (1994)  and  Yang  and  Liou  (1995a),  and  the  references  cited  therein. 
Significant  research  on  solving  light  scattering  by  large  hexagonal  ice  crystals  has  been  carried  out  in  the  past 
two  decades  by  means  of  the  geometric  ray-tracing  technique,  commonly  employed  to  identify  the  optical 
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phenomena  occurring  in  the  atmosphere.  However,  the  laws  of  geometric  optics  are  applicable  to  the  scattering 
of  light  by  an  ice  crystal  if  its  size  is  much  larger  than  the  incident  wavelength. 

In  recent  years,  we  have  developed  a  Monte  Carlo/geometric  ray-tracing  method  for  the  computation 
of  the  scattering,  absorption,  and  polarization  properties  of  ice  crystals  with  various  regular  and  irregular 
structures,  including  solid  and  hollow  columns,  single  and  double  plates,  dendrites,  bullet  rosettes,  and 
aggregates  (Takano  and  Liou  1995).  The  shapes  of  these  ice  crystals  are  defined  by  appropriate  geometric 
models  and  incident  coordinate  systems.  The  incident  photons  are  traced  with  a  hit-and-miss  Monte  Carlo 
method  and  followed  by  geometric  reflection  and  refraction  on  the  crystal  boundaries  in  which  absorption  can 
be  accounted  for  by  means  of  stochastic  procedures.  Phase  functions  for  shapes  covering  all  possible  ice 
crystal  structures  can  be  computed  using  this  method.  Thus,  at  this  point,  we  can  compute  the  single-scattering 
properties  for  large  ice  crystals  of  any  shape  that  can  be  defined  mathematically  or  numerically. 

The  geometric  ray-tracing  requires  the  principle  of  localization  in  which  the  particle  size  must  be  larger 
than  the  incident  wavelength.  In  addition,  it  is  assumed  that  the  energy  attenuated  by  the  scatterer  may  be 
decomposed  into  equal  extinction  from  diffraction  and  Fresnelian  rays  so  that  the  extinction  efficiency  is  equal 
to  2  regardless  of  the  particle  size  parameter.  To  circumvent  a  number  of  shortcomings  in  the  geometric  ray- 
tracing  approach,  we  have  developed  a  novel  improvement  by  mapping  the  equivalent  tangential  electric  and 
magnetic  currents  on  the  particle  surface  that  are  determined  from  geometric  reflection  and  refraction  to  fer- 
field  by  means  of  the  fundamental  electromagnetic  wave  theory  (Yang  and  Liou  1996a). 

For  verification  of  the  limitation  of  the  improved  geometric  ray-tracing  method,  referred  to  as  the 
geometric-optics/integral-equation  method  (GOM2),  we  have  also  developed  a  finite-difference  time  domain 
(FDTD)  technique,  pioneered  by  the  electrical  engineers  for  the  identification  of  irregular  objects.  It  is  a 
numerical  technique  for  the  solution  of  the  Maxwell  equations  using  appropriate  absorbing  boundary  conditions 
and  is  considered  to  be  the  "exact"  numerical  solution  for  light  scattering  by  particles,  as  verified  by  the  exact 
Mie-type  results  for  long  circular  cylinders  and  spheres  (Yang  and  Liou  1996b).  Because  of  solution 
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convergence  and  required  computer  time,  the  FDTD  method  is  applicable  to  size  parameters  smaller  than  about 
20.  The  conventional  geometric  optics  method  breaks  down  when  the  size  parameters  of  ice  crystals  are 
smaller  than  about  40-100,  depending  on  whether  the  computations  are  for  the  cross  sections  or  phase  matrix. 
Based  on  comparison  with  the  results  computed  from  the  FDTD  method,  the  GOM2  technique  is  shown  to  be 
applicable  to  ice  crystal  size  parameters  on  the  order  of  15-20.  By  combining  the  GOM2  and  FDTD  methods, 
referred  to  as  the  unified  theory  for  light  scattering,  calculations  for  the  scattering  and  absorption  properties 
of  ice  crystals  for  all  sizes  and  shapes  can  be  carried  out. 

Figures  3-5  show  the  computed  phase  functions  for  the  three  representative  size  distributions  described 
in  Section  2.2  and  for  the  75  ARES  channels.  Their  associated  single-scattering  parameters  are  given  in  Tables 
1-3.  In  the  calculations,  it  is  assumed  that  the  ice  crystal  size  distributions  are  composed  of  50%  bullet 
rosettes,  30%  hollow  columns  and  20%  plates.  It  is  noted  that  between  2  to  2.5  pm,  the  22°  and  46°  halos, 
and  the  backscattering  peak  are  prominent  for  the  December  5  case,  but  are  less  distinct  for  the  other  two  cases. 
For  the  cold  cirrus,  significant  side  scattering  is  evident.  In  addition,  for  the  cold  cirrus,  the  forward  scattering 
peaks  are  one  to  two  orders  of  magnitude  smaller  than  those  for  the  other  two  size  distributions. 

2.4  Effects  of  Small  Ice  Crystals  on  the  Single-Scattering  Properties 

Recent  observations  reveal  that  there  may  be  a  large  amount  of  small  ice  crystals  on  the  order  of  10 
pm  in  cirrus  clouds  (Platt  et  al.  1989).  Small  ice  crystals  were  also  observed  in  abundant  amount  in  tropical 
convective  anvils  (McFarquhar  and  Heymsfield  1996;  Heymsfield  and  McFarquhar  1996)  as  well  as  in  contrail 
cirrus  (Gayet  et  al.  1996)  and  in  orographic  wave  clouds  (Heymsfield  and  Miloshevich  1993).  The  presence 
of  small  ice  crystals  has  been  a  subject  of  constant  debate  in  the  study  of  cirrus  microphysical  and  radiative 
properties.  Due  to  the  high  scattering  efficiency  of  small  ice  crystal,  the  reflectance  of  cirrus  clouds  is  very 
sensitive  to  them.  To  accurately  quantify  the  effect  of  small  ice  crystals  on  the  cirrus  cloud  radiative 
properties,  reliable  ice  crystal  size  distributions  for  the  size  range  between  ~5  and  100  pm  are  needed.  Until 
recently,  it  was  extremely  difficult  to  measure  the  sizes  of  small  ice  crystals  with  confidence.  The  FSSP  that 
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Figure  3  Phase  functions  for  the  December  5  ice  crystal  size  distribution 
and  for  the  ARES  75-channel  spectral  wavelengths. 
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Figure  3  Continued. 


-  2.575|jm 

.  2.625pm 

-  2.675pm 

—  -  •  2.725pm 


Scattering  Angle  ( 0 ) 


0  60  120 
Scattering  Angle  ( ° ) 


Figure  4  Same  as  Figure  3,  except  for  cold  cirrus. 
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Figure  5  Same  as  Figure  3,  except  for  warm  cirrus  clouds. 
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Figure  5  Continued. 
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Figure  5  Continued. 
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Table  1 --continued. 
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Table  2.Same  as  Table  1,  except  for  the  cold  cirrus  ice  crystal  size  distribution. 
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Table  2-continued. 
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Table  3.  Same  as  Table  1,  except  for  the  warm  cirrus  ice  crystal  size  distribution. 
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Table  3--continued. 
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has  been  used  for  quantifying  the  spherical  droplet  size  distribution  has  not  been  calibrated  for  deriving  the 
nonspherical  ice  crystal  size  distributions  between  ~5  and  100  pm.  On  the  other  hand,  a  conservative  estimate 
of  the  lower  limit  of  ice  crystal  size  distributions  from  2D-C  probes  is  around  75  pm.  Given  that  there  are 
uncertainties  in  the  ice  crystal  size  distributions  obtained  from  either  FSSP  or  2D-C  probe,  two  fundamental 
aspects  need  to  be  explored.  First,  do  we  need  to  incorporate  ice  crystal  size  distribution  for  sizes  less  than  75 
pm  in  the  light  scattering  calculations?  If  the  answer  is  yes,  then  what  would  be  the  optimal  cut-off  size  for 
the  small  ice  crystal  size  distribution?  Second,  what  is  the  sensitivity  of  the  single-scattering  parameters  to  the 
uncertainties  in  the  small  ice  crystal  size  distribution? 

To  address  the  preceding  two  questions,  we  have  used  the  size  distributions  obtained  by  FSSP,  1D-C 
and  2D-C  probes  on  board  the  HARP  during  the  Hanscom  AFB  mission  flight  on  September  16,  1995,  for  light 
scattering  calculations.  The  FSSP  sizes  the  ice  crystals  from  1  to  96  pm;  the  1D-C  probe  samples  ice  crystals 
from  20  to  300  pm;  and  the  2D-C  probe  measures  the  maximum  dimensions  of  the  projection  of  ice  crystals 
from  25  to  800  pm.  Examination  of  these  ice  crystal  size  distributions  reveals  that  the  size  distributions 
derived  from  the  1D-C  and  2D-C  probes  basically  agree  with  each  other  in  terms  of  the  functional  shapes  and 
size  ranges  for  the  size  distribution,  although  the  differences  in  some  size  bin  concentrations  can  be  as  large 
as  one  order  of  magnitude.  Thus,  we  decided  to  use  the  2D-C  derived  ice  crystal  size  distribution  for  our 
purpose  because  its  size  range  covers  typical  ice  crystal  size  ranges.  Moreover,  we  note  that  the  FSSP  derived 
size  distributions  appear  to  be  unreliable,  because  the  size  bin  concentrations  are  unrealistically  high  for  sizes 
less  than  20  pm.  Thus,  to  obtain  the  ice  crystal  size  distribution  for  sizes  less  than  75  pm,  we  use  the  Gamma 
function  to  extrapolate  the  2D-C  size  distribution  to  1  pm. 

In  this  study,  three  composite  ice  crystal  size  distributions  obtained  by  the  2D-C  probe  are  constructed. 
They  are  named  as  cases  1  to  3,  and  are  plotted  in  Figure  6.  The  portion  of  the  ice  crystal  size  distribution 
larger  than  75  pm  was  derived  from  the  2D-C  probe  measurements.  The  size  distributions  between  1  and  75 
pm  were  obtained  from  the  extrapolation.  The  upper  bound  of  the  maximum  dimension  is  about  600  pm. 
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Figure  6  The  ice  crystal  size  distributions  based  on  HARP  data  (09/16/95)  with  the 
lower  limit  of  75|im  and  extrapolations  to  20pm  and  1pm  based  on  the 
Gamma  function  and  an  optimization  search  method. 
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There  was  no  specific  information  on  the  fractional  distribution  of  particular  ice  crystal  shapes. 

To  study  the  radiative  effects  of  small  ice  crystals,  we  use  the  unified  theory  for  light  scattering  by  ice 
crystals  to  compute  the  single-scattering  properties  of  the  aforementioned  three  ice  crystal  size  distributions, 
including  the  phase  function  (asymmetry  factor),  single-scattering  albedo,  and  extinction  coefficient.  To 
investigate  the  effects  of  different  cut-off  sizes  on  the  single-scattering  properties,  three  cut-off  sizes  are 
prescribed:  1  pm  20  pm,  and  75  pm.  The  cut-off  size  at  75  pm  corresponds  to  no  extrapolation  at  all.  Table 
4  lists  the  single-scattering  parameters  at  2.8  pm  for  the  three  cut-off  sizes  for  each  of  the  ice  crystal  size 
distributions.  Figure  7  shows  the  associated  normalized  phase  function  for  these  three  cases.  It  is  assumed 
that  ice  crystals  are  composed  of  randomly  oriented  ice  columns  and  plates.  The  extinction  coefficient  is 
significantly  affected  by  the  variation  of  cut-off  sizes.  The  value  for  cut-off  at  1  pm  is  virtually  the  same  as 
that  for  cut-off  at  20  pm  for  all  three  cases,  but  the  extinction  coefficient  for  cut-off  at  75  pm  is  much  smaller 
than  that  for  cut-off  at  20  pm.  Similar  patterns  are  displayed  for  the  associated  normalized  phase  functions. 
The  normalized  phase  functions  for  cut-off  at  1  pm  is  virtually  the  same  as  those  for  cut-off  at  20  pm  for  all 
three  cases,  but  the  normalized  phase  functions  for  cut-off  at  75  pm  is  much  smaller  than  those  for  cut-off  at 
20  pm.  Thus,  it  is  important  to  take  into  account  the  small  crystal  size  distribution  between  20  and  75  pm  in 
the  light  scattering  calculations. 

We  also  examined  the  effects  of  variation  in  the  small  ice  crystal  number  concentration  on  the  single¬ 
scattering  properties.  We  use  the  case  1  size  distribution  for  this  purpose.  Figure  8  shows  the  log-log  linearly 
extrapolated  ice  crystal  size  distribution.  A  log-log  linear  extrapolation  was  applied  to  the  distribution  based 
on  the  FSSP  data.  The  symbol  x  denotes  the  number  concentration  corresponding  to  the  mean  ice  crystal  size 
(15  pm)  determined  from  the  FSSP  data.  The  symbols  y  and  z  are  the  number  concentrations  at  twice  and  half 
the  value  of  x,  respectively.  The  symbol  s  represents  the  number  concentration  for  the  particle  size  of  75  pm. 
Three  extrapolation  curves  are  subsequently  obtained  as  the  straight  line  segments  linking  s  and  x  (mean),  s 
and  y  (mean*2),  and  s  and  z  (mean/2). 
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Table  4.  The  single-scattering  parameters  for  the  ice  crystal  size  distributions  depicted  in  Fig.6  at  X=2.8|im. 
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ire  7  The  normalized  phase  functions  for  the  ice  crystal  size  distributions  given  in  Figure  6 
at  X=2.8pm,  computed  from  a  modified  ray-tracing  method  for  randomly  oriented 
hexagonal  columns  and  plates. 
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Figure  8  The  ice  crystal  size  distributions  based  on  2D-C  and  FSSP  measurements 
during  HARP  mission  flight  on  September  16,  1995  (Case  1).  The  lower 
limit  of  the  2D-C  size  distribution  is  75  pm.  Three  log-log  linear 
extrapolations  from  75pm  to  20pm  are  displayed  ( see  text). 


Table  5  lists  the  single-scattering  parameters  at  2.8  pm  for  the  ice  crystal  size  distributions  obtained 
from  the  extrapolated  values.  Figure  9  shows  the  associated  normalized  phase  function  for  the  three 
extrapolated  curves.  Similar  to  the  study  of  cut-off  sizes,  only  the  extinction  coefficient  is  affected  by  the 
variation  of  small  ice  crystals.  The  single-scattering  albedo,  asymmetry  factor,  and  the  phase  function  do  not 
vary  significantly,  however.  From  the  preceding  studies,  we  conclude  that,  to  obtain  a  correct  value  for  the 
extinction  coefficient,  it  is  necessary  to  include  a  reliable  extrapolated  ice  crystal  size  distribution  between  20 
and  75  pm  in  the  light  scattering  calculations. 

3.  RETRIEVAL  OF  CIRRUS  CLOUD  PARAMETERS  USING  ARES  5.1-5.3  AND  3.7  pm 
CHANNEL  DATA 

In  this  section,  we  present  an  airborne  retrieval  algorithm  to  infer  cirrus  cloud  temperature,  optical 
depth,  and  ice  crystal  mean  effective  size  using  the  ARES  IR  spectrometer  data.  The  retrieval  method  utilizes 
the  5. 1-5. 3  pm  water  vapor  band-wing  radiances  to  compute  cirrus  cloud  temperatures  and  IR  emissivities, 
from  which  the  associated  optical  depth  is  derived  based  on  radiative  transfer  parameterizations.  The  ice 
crystal  mean  effective  size  is  obtained  by  matching  the  computed  sum  of  the  3.7  pm  solar  reflected  and  IR 
emitted  radiances  with  those  from  the  ARES  Ch.  42  (3.755  pm)  radiances.  The  selection  of  these  channels  is 
based  on  the  examination  of  the  computed  weighting  functions  of  a  midlatitude  summer  atmosphere  and  the 
inspection  of  images  of  ARES  channel  radiances  over  the  WB-57F  flight  tracks.  We  applied  the  retrieval 
algorithm  to  the  ARES  data  collected  during  a  mission  flight,  which  was  carried  out  on  board  the  WB-57F  on 
September  16, 1995  over  the  western  Boston  area.  The  retrieved  cloud  parameters  are  compared  with  the  in- 
situ  2D-probe  measurements  made  from  the  Aeromet's  Learjet  as  well  as  cloud  heights  derived  from  ground- 
based  8.6-mm  radar  returns. 

Cloudy  radiances  have  been  measured  by  the  IR  imaging  spectrometer/radiometer  on  board  the  ARES 
Programs'  WB-57F  high-altitude  research  aircraft  operated  by  the  NASA  Johnson  Space  Flight  Center.  The 
WB-57F  is  capable  of  reaching  altitudes  in  excess  of  60,000  ft,  and  carries  a  crew  of  two,  permitting 
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Table  5 .  The  single-scattering  parameters  for  the  ice  crystal  size  distributions  depicted  in  Figure  8  at 

A  =  2.8  pm. 


mean  FSSP/2 

mean  FSSP 

mean  FSSP  *  2 

Extinction 
Coefficient  (1/km) 

0.407890 

0.446238 

0.522934 

Single¬ 
scattering  Albedo 

0.526002 

0.528148 

0.531495 

Asymmetry  factor 

0.985027 

0.984543 

0.983796 
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Normalized  Phase  Function 


Figure  9  The  normalized  phase  function  for  the  randomly  oriented  ice 

columns  and  plates  at  X=2.8pm.  The  ice  crystal  size 
distributions  of  ice  crystals  are  depicted  in  Figure  8. 
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instrument  operation  by  a  dedicated  flight  crew.  The  sensor  was  mounted  on  the  lower  fuselage  surface,  with 
downward  viewing  while  in  level  flight.  The  ARES  sensor  for  the  present  study  has  been  operated  as  a  75- 
channel  imaging  spectrometer.  It  is  operated  in  a  pushbroom  mode  with  45-pixel  array  pointing  in  the  nadir 
direction  taking  data  at  the  standard  rates  of  10,  20,  40,  or  80  scans/sec.  The  optical  elements  of  ARES 
provide  a  pixel  IFOV  of  slightly  more  than  1  mrad,  which  corresponds  to  about  15  m  ground  resolution  for  a 
flight  altitude  of  15  km.  The  spacing  of  the  ARES  spectral  wavelengths  appears  to  depend  on  the  optical 
arrangement  with  the  channel  band  widths  that  vary  between  0.02  and  0.07  pm. 

For  validation  of  the  retrieved  cloud  temperature/height,  we  have  obtained  the  cloud  boundary  height 
data  from  the  collocated  and  coincident  ground-based  8.6-mm  TPQ-1 1  radar  measurements.  For  comparison 
of  the  retrieved  cloud  optical  depth  and  mean  effective  sizes  with  in-situ  measurements,  we  have  also  acquired 
the  meteorological  and  ice  microphysical  data  collected  from  the  High  Altitude  Reconnaissance  Platform 
(HARP),  which  is  a  Learjet  36A,  flying  in  conjunction  with  the  ground-based  radar  and  ARES  measurements. 
The  HARP  carried  instruments  for  taking  both  in-situ  and  remote  atmospheric  measurements.  The  in-situ 
instruments  included  FSSP,  1D-C  and  2D-C  probes  for  particle  size  distribution  measurements.  Ambient 
temperatures  were  measured  by  the  Rosemount  Total  Temperature  Sensor  with  an  uncertainty  of  ±1  K. 

3.1  Principle  of  Retrieving  Cirrus  Temperature,  Emissivity,  and  Ice  Crystal  Size  Using  ARES  Channels 

The  present  retrieval  method  for  inferring  cirrus  temperature  and  emissivity  is  based  on  the  numerical 
solution  of  the  following  set  of  algebraic  equations  (Liou  et  al.  1990;  Ou  et  al.  1993,  1995)  as  follows: 

Ri  =(l-ei)Rai  +  €i  B;(TC),  (1) 

e;  =  1  -  exp(-kjT  ),  (2) 

where  R,  is  the  upwelling  radiances  at  the  top  of  the  cirrus  cloudy  atmosphere  for  the  ith  selected  ARES 
channel;  R^  is  the  corresponding  radiances  reaching  the  cloud  base  for  the  ith  channel;  e;  is  the  cloud  IR 
emissivities,  Tc  is  the  mean  cirrus  temperature;  B;  (Tc)  is  the  Planck  intensities  at  Tc,  t  is  the  visible  optical 
depth  and  lq  is  the  effective  extinction  coefficient  that  has  been  determined  from  the  calculated  single-scattering 
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properties  based  on  a  number  of  measured  ice  crystal  size  distributions.  The  behavior  of  the  cloud  emissivity 
in  terms  of  its  dependence  on  wavelength  and  cloud  parameters  has  been  examined  previously  (Liou  et  al. 
1990). 

By  straightforward  algebraic  manipulations,  Eqs.  (1)  and  (2)  for  a  pair  of  channels  (i=l  and  2)  can 
be  combined  into  a  single  equation  in  the  form 

'VWl  .  [R,  -B.dJlMc, 

Rtf-JMTjJ  [R„  -  B,  (T*)] 

IfR,,  and  are  known,  B2(TC)  and  B,(TC)  are  linearly  correlated,  and  k2/ki  ~  1.0,  Eq.  (3)  becomes  a  linear 
equation.  Numerical  solution  of  linear  equations  is  relatively  straightforward  and  easy  to  program,  as 
compared  to  the  solution  of  nonlinear  Eqs.  (l)-(2). 

The  effective  extinction  coefficients  for  the  retrieval  channels  have  been  computed  using  die 
adding/doubling  program  in  conjunction  with  the  development  of  this  remote  sensing  algorithm.  Eleven  ice 
crystal  size  distributions  from  aircraft  in-situ  measurements  have  been  chosen  (Fu  and  Liou  1993).  Their  mean 
effective  ice  crystal  sizes  range  between  23.9  and  123.6  pm.  The  single-scattering  properties  including  phase 
functions,  single  scattering  albedos,  and  extinction  coefficients  have  been  computed  for  each  size  distribution 
and  for  each  of  the  selected  ARES  channels  using  the  unified  theory  for  light  scattering  by  ice  crystals.  To 
optimize  the  computational  effort,  the  forward-scattering  peak  of  each  phase  function  was  truncated  following 
the  procedure  developed  by  Takano  and  Liou  (1989).  Radiances  for  each  channel  and  for  a  range  of  optical 
depths  and  ice  crystal  sizes  were  computed  using  the  aforementioned  adding/doubling  program.  The  cirrus 
cloud  emissivity,  which  is  defined  as  the  ratio  of  the  the  upwelling  radiance  to  the  cloud  Planck  function,  was 
obtained  as  functions  of  the  optical  depth,  mean  effective  size,  and  channel  wavelengths.  Finally,  the  effective 
extinction  coefficients  for  each  channel  and  for  each  size  distribution  was  obtained  as  the  slope  of  the  linear 
fitting  of  -ln(l-6i)  vs.  x.  Figure  10  shows  the  effective  extinction  coefficient  (k-value)  as  a  function  of  the 
mean  effective  size  for  eight  wavelengths.  For  all  wavelengths,  the  k-values  increase  with  the  mean  effective 
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Figure  10  The  k-value  as  a  function  of  De  for  the  retrieval  channels.  Eleven 
size  distributions  have  been  used  for  single-scattering 
computations.  The  single-scattering  properties  are  used  in  an 
adding/doubling  code  to  obtain  the  cloud  emissivity. 


size.  This  functional  dependence  of  the  k-values  has  been  incorporated  into  the  retrieval  program.  More 
importantly,  the  diagram  shows  that  the  values  of  the  three  pairs  of  k-ratios  —  k(5 . 1 89  pm)/k(5 .116  pm), 
k(5.261  pm)/k(5.116  pm),  and  k(5.332  pm)/k(5.116  pm)— are  all  close  to  1.0  for  the  range  of  the  mean 
effective  size  considered  here.  This  is  one  of  the  requirements  for  the  linear  correlation  of  pairs  of  ARES 
channel  radiances. 

In  order  to  infer  the  cirrus  emissivity,  ei5  and  temperature,  Tc,  from  multiple  radiance  measurements 
in  cloudy  conditions,  the  upwelling  radiances  at  the  cloud  base,  R^,  must  be  known.  In  principle,  it  can  be 
generated  from  forward  radiation  transfer  calculations  if  the  surface  temperature,  temperature  and  humidity 
profiles,  and  cloud  base  height  are  given.  However,  a  more  direct  and  reliable  method  to  determine  R^  would 
be  to  find  clear  column  radiances  from  the  available  data. 

Assuming  that  the  water  vapor  absorption/emission  effects  are  negligible  above  cirrus  clouds,  we  may 
approximate  the  cloud-base  upwelling  radiance  using  the  measured  clear  radiance  at  the  ARES  flight  level. 
This  approximation  has  been  applied  to  the  retrieval  of  cirrus  cloud  parameters  using  AVHRR  3.7  and  10.9 
pm  window  bands  (Ou  et  al.  1993)  and  to  the  airborne  cirrus  retrieval  using  6.5  and  10.5  pm  radiometer  data 
(Liou  et  al.  1990).  Once  is  determined,  Eqs.  (l)-(2)  can  be  solved  by  a  numerical  iteration  method.  Since 
there  are  two  independent  parameters,  N-l  pairs  of  the  N  (>  1)  channel  data  are  used  in  the  retrieval  program. 
N-l  sets  of  the  retrieved  temperature  and  emissivity  are  then  averaged  to  obtain  their  respective  mean  values. 

Upon  retrieval  of  the  visible  optical  depth  and  the  cloud  temperature  using  thermal  IR  channels,  we 
can  obtain  the  mean  effective  size  using  the  observed  3.7  pm  channel  radiances.  We  first  assume  a  trial  ice 
crystal  size.  The  3.7  pm  thermal  IR  radiance  can  be  calculated  according  to  Eq.  (1).  The  cloud-base  upwelling 
3.7  pm  IR  radiance  can  be  obtained  by  converting  the  4.75  pm  clear-pixel  radiance  into  its  equivalent  Planck 
function  at  3.7  pm,  since  both  channels  are  in  the  window  spectral  region.  The  3.7  pm  solar  reflectance  is 
determined  from  look-up  tables,  which  are  generated  for  the  same  six  size  distributions  as  described  in  Rao  et 
al.  (1995)  and  for  appropriate  solar-viewing  geometries,  using  the  adding/doubling  method  developed  by 
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Takano  and  Liou  (1989).  Finally,  the  total  3.7  pm  radiance  is  calculated  as  the  sum  of  the  thermal  IR  and 
solar  radiances.  The  radiance  values  thus  obtained  are  subsequently  compared  with  the  observed  3.7  pm 
radiances.  The  effective  mean  ice  crystal  size  can  then  be  obtained  iteratively  by  adjusting  the  trial  value  based 
on  the  difference  between  the  computed  and  observed  radiances.  Pixels  having  a  retrieved  mean  effective  size 
greater  than  123.6  pm  or  less  than  23.9  pm  are  discarded.  Figure  1 1  gives  a  schematic  flow  diagram  for  the 
retrieval  of  cirrus  cloud  temperature,  emissivity  (optical  depth)  and  mean  ice  crystal  size. 

3.2  Application  of  the  Retrieval  Scheme  to  ARES  Data 

The  preceding  cirrus  cloud  retrieval  scheme  is  applied  to  the  ARES  spectrometer  data  obtained  during 
a  WB-57F  flight  that  was  launched  on  September  16,  1995.  Data  collected  for  Tracks  1,  3,  5,  and  8  of  that 
mission  were  analyzed.  For  Track  1,  the  WB-57F  flew  between  2055  and  2100  at  a  cruising  altitude  of  15 
km  along  the  southeast-northwest  direction  (Figure  12).  The  WB-57  flew  directly  over  Hanscom  Air  Force 
Base  (represented  by  the  symbol  "x"),  where  the  ground-based  8.6-mm  radar  was  located.  During  the  period 
of  Tracks  1  and  3,  HARP  was  also  flying  below  the  WB-57F,  although  the  two  aircraft  were  not  exactly 
collocated. 

3.2.1  ARES  Channel  Spectral  Properties 

Adequate  implementation  of  the  retrieval  scheme  requires  knowledge  of  the  ARES  channel  spectral 
properties  and  the  adequate  selection  of  useful  ARES  channel  data.  Figure  13  presents  the  upwelling 
radiances  associated  with  solar  reflection  and  thermal  emission  in  the  2-6.4  pm  region  to  assess  the  relative 
importance  of  the  two  radiation  sources.  For  the  solar  reflection,  we  have  used  the  updated  solar  xrradiance 
data  provided  for  our  use  by  G.  Anderson.  The  results  are  shown  for  the  cosine  of  two  solar  zenith  angles  and 
for  a  surface  reflection  of  0.1.  For  the  latter,  Planck  function  curves  for  three  temperatures  covering  the  range 
of  relevant  earth's  atmospheres  are  shown.  It  is  clear  that  for  3  pm,  the  solar  reflection  is  the  primary 
contributor  to  radiation,  while  for  4.5  pm,  the  thermal  IR  emission  is  the  predominant  radiation  source. 
Between  these  two  wavelengths,  the  relative  importance  of  the  solar  reflection  and  thermal  IR  emission  is 
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Figure  13  Upwelling  radiances  for  solar  reflection  and  thermal  emission  in  the  2-6.5  pm  region.  Fx  denotes  the  solar  irradiances, 
|iD  is  the  consine  of  solar  zenith  angle,  a  is  the  surface  albedo,  and  Bx  is  the  Planck  function. 


dependent  on  the  sun's  position,  the  surface  albedo,  and  the  atmospheric  thermodynamic  state.  For  the  radiative 
transfer  simulation  and  the  remote  sensing  algorithm  development,  both  components  must  be  accounted  for  in 
this  intermediate  spectral  range. 

In  the  2-6.4  pm  spectral  region,  several  absorption  bands  exist.  These  include  the  2.7  pm  H20  band 
produced  by  v,  and  v2  fundamentals,  the  3.2  pm  H20  band  associated  with  2v2,  the  4.3  pm  C02  due  to  v3 
fundamental  and  P  and  R  branches,  and  the  H20  6.3  pm  vibrational-rotational  band.  There  are  three  windows 
located  at  the  2.2,  3.7,  and  4.7  pm  regions.  The  spacing  of  the  ARES  spectral  wavelengths  (wavenumbers) 
appears  to  depend  on  the  optical  arrangement,  but  75  channels  are  available  for  cloud  detection.  Moreover, 
solar  reflection  contributes  significantly  to  radiances  for  wavelengths  shorter  than  5  pm.  Figure  14  shows  a 
typical  example  of  the  ARES  75-channel  spectra  for  clear  and  cloudy  areas.  Respective  spectral  regions  for 
the  absorption  bands  are  noted  in  the  diagram.  Large  measurement  noises  exist  for  Chs.  1-25  (1.9-2.7  pm), 
so  that  radiances  for  these  channels  are  not  useful.  For  Chs.  26-75,  the  aforementioned  four  absorption  bands 
and  two  window  bands  are  clearly  displayed  as  local  minima  and  maxima. 

3.2.2  Selection  of  Retrieval  Channels 

To  select  a  number  of  appropriate  channels  for  the  cirrus  cloud  retrieval,  we  have  examined  all  the 
ARES  channel  weighting  functions  and  images.  Typical  weighting  functions  were  generated  using  the  vertical 
transmittance  profiles  obtained  from  MODTRAN  calculations  employing  a  midlatitude  summer  temperature 
and  humidity  profiles.  Figure  15  shows  the  weighting  functions  for  all  the  75  ARES  channels.  There  are 
features  that  correspond  to  the  distribution  of  absorption  bands.  Over  the  window  band  regions  (2.2  pm,  3.75 
pm,  and  4.7  pm),  the  weighting  function  peaks  at  or  near  the  surface.  Over  the  band  center  of  the  2.7  pm  H20, 
4.3  pm  C02,  and  6.3  pm  H20  absorption  bands,  the  weighting  function  peaks  in  the  upper  atmosphere.  Over 
the  band  wing  regions,  the  weighting  function  peaks  at  intermediate  levels.  These  physical  properties  are  the 
basis  for  the  selection  of  retrieval  channels.  In  principle,  it  is  most  ideal  to  choose  channels  with  the  weighting 
functions  peaking  in  the  lower  atmosphere.  For  channels  with  the  weighting  functions  peaking  in  the  upper 
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Figure  14Typical  example  of  the  ARES  75-channel  spectra  for  clear  and 
cloudy  areas. 
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atmosphere,  the  radiances  tend  to  be  dominated  by  gaseous  emission  from  the  upper  atmosphere,  thus 
suppressing  the  cloud  signals.  On  the  other  hand,  for  channels  with  the  weighting  functions  peaking  at  the 
surface,  inhomogeneity  can  cause  a  large  spread  in  clear  radiances  and  thus  lead  to  inaccurate  cloud  retrieval 
results. 

Moreover,  we  can  also  use  radiance  images  to  facilitate  the  selection  of  channels.  We  have  acquired 
from  SC1TEC  two  CDROMs  containing  the  ARES  data  16-9-95-1A,  the  quick-look  software  ADAPT  from 
SC1TEC  which  allows  us  to  display  graphically  the  large  volume  of  the  ARES-SSDS  dataset  on  SGI,  and  the 
decoding  software  "hdc2bip"  which  enables  us  to  decode  the  compressed  data  set.  Using  ADAPT,  we  were 
able  to  scan  the  image  of  the  entire  flight  track  for  each  ARES  channel,  and  to  select  a  scene  which  displays 
a  distinct  contrast  between  high  (clear)  and  low  (cloudy)  radiance  areas.  Figure  16  shows  the  radiance  images 
for  Chs.  41(3.681  pm),  42(3.755  pm),  55(4.747  pm),  60(5.116  pm),  and  63(5.332  pm)  over  the  selected 
scene.  The  first  three  channels  are  within  the  3.75  pm  and  4.7  pm  window  regions,  while  the  last  two  channels 
are  in  the  6.3  pm  water  vapor  band-wing  region.  For  each  image,  the  horizontal  and  vertical  scales  are  -500 
m  and  -2.5  km,  respectively. 

From  these  images,  it  is  noted  that  the  upper  1/3  of  the  domain  is  clear,  while  the  rest  of  the  domain 
is  cloudy.  Over  the  clear  area,  the  images  for  the  three  window  channels  show  significant  variations  of  the 
emitted/  reflected  radiances  due  to  the  inhomogeneity  of  the  surface  temperature  and  emissivity.  However,  for 
the  two  band-wing  channels,  radiances  over  the  clear  area  appear  to  be  relatively  uniform,  as  compared  to  the 
clear  radiances  for  die  window  channels,  because  the  significant  water  vapor  absorption/emission  effect  in  the 
lower  troposphere  blocks  the  horizontally  inhomogeneous  surface  emission.  The  narrow  distribution  of  the 
clear  radiances  is  a  desirable  feature  for  application  of  the  remote  sensing  algorithm.  A  smaller  spread  in  clear 
radiances  would  lead  to  less  errors  in  prescribing  the  mean  clear  radiances,  thus  increasing  the  retrieval 
accuracy  (Rao  et  al.  1995). 

Over  the  cloudy  area,  it  is  evident  that  the  image  for  the  4.7  pm  channel,  as  well  as  the  images  for  the 
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Figure  16  Radiance  images  for  ARES  channels  41,  42,  55,  60,  and  63  over  a  selected 
scene  with  strong  cloud/clear  contrast  within  the  Track  1  of  September  16 
flight. 


two  band-wing  channels,  contain  distinct  cloud  signatures.  For  the  two  3.7  pm  window  channels,  there  is  no 
distinguishable  cloud  signature  due  to  the  compensating  effects  of  strong  reflection  and  weak  emission  by  cloud 
particles.  Radiance  images  for  channels  at  the  absorption  band  center  (not  shown  here)  exhibit  no  cloud 
signature.  Based  on  the  images  shown  in  Figure  16,  together  with  the  weighting  function  profiles  shown  in 
Figure  15,  we  select  Ch.  60,  61  (5.189  pm),  62  (5.261  pm)  and  63  for  the  development  of  a  cirrus  cloud 
retrieval  algorithm.  Images  for  these  channels  contain  definitive  cloud  signatures  and  display  narrow  data 
spread  of  clear  radiances. 

3.2.3  Determination  of  Mean  Clear  Radiances 

To  retrieve  the  cloud  temperature  and  emissivity,  the  mean  clear  radiance  for  each  channel  within  the 
selected  domain  must  be  known.  Thus,  a  clear/cloud  detection  scheme  was  developed  to  identify  clear  data 
points.  It  is  based  on  the  physical  properties  of  the  ARES  equivalent  brightness  temperatures  for  selected 
channels.  In  this  scheme,  clear  pixels  are  separated  from  the  cloudy  pixels  by  applying  the  threshold  tests  to 
the  equivalent  brightness  temperature  data  for  all  selected  channels.  If  the  equivalent  brightness  temperatures 
for  given  pixels  are  greater  than  the  prescribed  threshold  values  for  all  selected  channels,  then  those  pixels  are 
categorized  as  clear.  Otherwise,  the  pixel  is  flagged  as  cloudy.  A  schematic  flow  diagram  depicting  the 
clear/cloudy  detection  scheme  is  given  in  Figure  17.  The  number  of  cloudy  pixels  may  be  overestimated  in  this 
scheme.  However,  by  doing  so,  we  can  ensure  that  the  detected  clear  pixels  are  truly  clear  pixels,  and  that  any 
cloudy  pixels  will  not  be  missed  when  the  retrieval  is  carried  out. 

The  mean  clear  radiance  for  each  channel  is  determined  based  on  the  average  of  radiances  of  all 
detected  clear  pixels  within  the  selected  scene  (- 10  km  strip).  To  determine  the  appropriate  threshold  values 
for  all  selected  rhannds  we  perform  statistical  analyses  on  the  radiance  data  over  the  domain  of  Figure  16  that 
has  distinct  clear-cloud  contrast.  Figures  1 8(a)-(d)  show  the  histograms  of  the  brightness  temperatures  for  Chs. 
60-63  within  the  domain  of  Figure  16.  Each  histogram  displays  a  distinct  bi-modal  shape.  The  high- and  low- 
radiance  modes  represent  the  mean  radiances  for  the  clear  and  cloudy  conditions,  respectively.  Because  these 
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Figure  17  Schematic  flow  diagram  depicting  the  clear/cloud  detection  scheme. 
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modes  are  so  distinct,  the  threshold  value  for  each  channel  can  be  easily  determined  as  the  radiance  that  is 
associated  with  the  minimum  frequency  between  the  two  modes.  Figures  19(a)-(e)  show  the  clear  (black)  and 
the  cloudy  (white)  pixels  detected  by  individual  and  combined  threshold  tests.  For  Chs.  60  and  61,  because 
of  their  rather  transparent  nature,  certain  pixels  with  higher  brightness  temperatures  are  identified  as  "clear". 
However,  these  "clear"  pixels  are  not  detected  by  the  other  two  channels.  For  this  reason,  the  combined 
threshold  tests  detect  more  cloudy  pixels  than  those  determined  from  individual  channels. 

3.2.4  Application  of  the  Retrieval  Scheme  to  ARES  Channel  Data 

Before  application  of  the  retrieval  scheme  to  the  ARES  channel  data,  it  would  be  interesting  to  examine 
whether  the  ARES  IR  channel  data  show  characteristics  of  linearity  in  the  two-channel  radiance  correlation 
[Eq.  (3)],  and  to  determine  whether  the  ARES  3.7  pm  channel  data  fall  within  the  range  of  the  computed 
values  so  that  these  data  can  be  used  for  the  retrieval  of  ice  crystal  mean  effective  sizes.  Figure  20  shows 
examples  of  the  scatter  diagram  for  three  pairs  of  ARES  channels  :  61-60,  62-60,  and  63-60.  In  each  frame, 
the  straight  line  denotes  the  correlation  of  Planck  functions  for  the  two  channels,  and  the  dots  are  radiance  data 
points  chosen  from  the  ARES  data  collected  during  the  September  16,  1995,  mission  flight  over  the  Hanscom 
Air  Force  Base  area,  where  a  strong  contrast  of  clear  and  cloudy  regions  was  present  at  the  flight  time  (see 
Section  3.2.2).  The  clustering  of  radiance  data  points  shows  a  distinctive  straight-line  pattern.  Basically,  there 
are  two  clusters  in  each  frame.  The  upper  cluster  corresponds  to  clear  data  points,  while  the  lower  longer 
clusters  are  associated  with  cloudy  radiances. 

Figure  21  shows  an  example  ofthe  retrieval  ofice  crystal  mean  effective  sizes  using  the  ARES  3.7  pm 
channel  data.  In  the  upper  panel,  the  computed  3.7  pm  solar  radiances  are  plotted  as  functions  of  optical  depth 
and  six  mean  sizes.  For  x  >  2,  these  radiances  are  not  sensitive  to  the  optical  depth,  but  are  quite  sensitive  to 
the  ice  crystal  size  distribution.  In  the  middle  panel,  the  computed  3.7  pm  IR  radiances  for  different  sizes  are 
plotted  against  the  optical  depth.  In  this  case,  the  cloud  temperature  is  250  K  and  the  clear  brightness 
temperature  is  291 K.  The  3.7  pm  IR  radiances  decrease  with  increasing  optical  depth  and  mean  effective  size. 
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Figure  20  Examples  of  scattere  diagram  for  three  pairs  of  ARES  channels:  60-61 , 60-62,  and  60-63. 
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Figure  21  Computed  3.7  pm  solar  radiances  (a),  IR  radiances  (b),  and  total 
radiances  (c)  as  functions  of  optical  depth  and  ice  crystal  size. 
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Finally,  in  the  lower  panel,  the  3.7  pm  total  radiances  are  plotted  against  the  optical  depth.  It  is  noted  that 
these  radiances  are  sensitive  to  both  the  optical  depth  and  ice  crystal  size.  We  display  the  ARES  measured  data 
collected  during  Track  3  of  the  September  16  flight  on  the  lower  diagram  and  show  that  they  lie  within  the  bulk 
range  of  the  computed  curves,  indicating  that  it  is  feasible  to  use  the  ARES  measured  data  to  retrieve  ice  mean 
effective  sizes. 

Applying  the  retrieval  scheme  to  the  detected  cloudy-pixel  radiance  data  of  the  selected  ARES 
channels,  we  can  determine  the  cirrus  cloud  temperatures,  emissivities  (optical  depth),  and  ice  crystal  sizes  for 
the  four  flight  tracks  of  September  16, 1995.  For  the  purpose  of  displaying  data,  we  obtain  the  array-averaged 
channel  equivalent  brightness  temperatures  and  the  retrieved  parameters  for  each  scan  within  the  flight  track. 
Variation  of  the  computed  array-averaged  Ch.  60  equivalent  brightness  temperatures,  cloud  temperatures,  and 
optical  depths  for  the  four  flight  tracks,  the  mean  effective  sizes  far  Tracks  1  and  3,  and  the  cloud  emissivities 
for  Tracks  5  and  8  along  the  flight  path  are  shown  in  Figures  22-25.  The  results  shown  in  these  figures  are 
based  on  the  application  of  the  retrieval  scheme  for  4  to  ~  10  scans  of  the  45-pixel  array  data.  However,  by 
using  an  efficient  numerical  scheme,  the  total  run  time  on  a  powerful  state-of-the-art  SGI  is  less  than  10  min 
for  each  track.  Retrieval  of  the  mean  effective  size  was  not  performed  for  Tracks  5  and  8,  because  these  tracks 
were  flown  after  nightfall  and  the  3.7  pm  radiances  are  dominated  by  the  thermal  IR  component,  which  is  less 
sensitive  to  the  ice  crystal  size  than  the  solar  component.  For  Track  1,  there  are  three  intervals  which  have 
been  identified  as  clear  (16-20  km,  25-28  km  and  30-32  km).  For  Track  3,  one  clear  period  is  identified 
between  16  and  20  km.  Tracks  5  and  8  show  no  clear  periods.  The  retrieved  cloud  temperatures  vary  between 
230  and  255  K.  The  retrieved  cloud  optical  depths  are  0-2,  0-4,  1.6-7,  and  1. 2-7.2  for  the  four  tracks.  The 
retrieved  mean  effective  ice  crystal  sizes  for  Tracks  1  and  3  are  generally  between  25  and  75  pm.  Based  on 
these  retrieval  results,  the  cirrus  under  study  was  relatively  warm  and  was  composed  of  small  ice  crystals.  The 
cirrus  clouds  for  Tracks  1  and  3  were  optically  thin,  but  they  became  thicker  during  the  last  two  flight  tracks. 
The  reliability  of  the  retrieved  results  is  partially  confirmed  by  the  following  validation  efforts. 
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Figure  22  Array  averaged  channel  equivalent  brightness  temperatures, 
retrieved  cloud  temperature,  cloud  optical  depth,  and  effective 
size  as  functions  of  flight  distance  for  Track  1 , 9/16/95. 


3.2.5  Validation  of  Retrieval  Results 

The  cirrus  clouds  on  September  16  was  also  remotely  observed  by  the  ground-based  8.6-mm  radar 
which  operated  continuously  during  all  four  flight  tracks.  Each  track  passed  the  radar  site  once.  Moreover, 
the  same  area  of  cirrus  clouds  were  also  sampled  in-situ  by  cloud  probes  on  board  HARP,  which  flew 
simultaneously  below  ARES  in  the  vicinity  of  the  ground-observation  site.  Because  of  the  coincidence  and 
collocation  of  the  ARES,  radar,  and  cloud  probes  operations,  we  can  validate  the  ARES  retrieved  cirrus  cloud 
temperatures/heights,  optical  depths  and  mean  effective  sizes  using  the  radar  and  cloud  probe  measurements. 
These  validation  efforts  are  described  in  the  following. 

To  validate  the  retrieved  cirrus  cloud  temperature/height,  we  use  the  time  series  of  cloud  boundary 
(cloud  top  and  base)  altitudes  derived  from  the  8.6-mm  radar  measurement,  together  with  the  upper 
tropospheric  temperature  measurement  by  HARP.  Figure  26(a)  shows  the  atmospheric  temperature  profile 
between  5  and  11  km.  These  temperature  values  were  obtained  between  2029  and  2341  UTC.  Within  the 
altitude  range,  the  temperature  varies  between  -5°  and  -55°C.  Figure  26(b)  shows  the  cloud  boundary 
altitudes  derived  from  the  8.6-mm  radar  measurements  between  2000  and  2400  UTC.  Within  this  time  period, 
the  8.6-mm  radar  observed  continuous  presence  of  a  single  layer  of  cirrus  cloud.  Between  2000  and  2100 
UTC,  the  observed  cloud  top  and  base  were  located  at  8  and  6  km,  respectively.  After  2100  UTC,  the  cloud- 
top  altitude  fluctuates  between  8  and  10  km,  while  the  cloud  base  altitudes  largely  remains  at  6  km,  except  for 
the  period  between  2130  and  2230  UTC.  We  convert  the  retrieved  cloud  temperatures  into  cloud  altitudes 
using  the  HARP-derived  sounding.  (There  was  no  surface  sounding  launched  around  the  flight  time  of  the  WB- 
57F.)  For  Tracks  1  and  3,  when  WB-57F  flew  over  the  radar  site,  the  area  happened  to  be  clear,  although  the 
radar  recorded  very  thin  cirrus.  Thus,  comparison  of  retrieved  cloud  altitudes  with  the  radar  data  was  not 
available.  For  Tracks  5  and  8,  the  retrieved  cloud  temperatures  over  the  radar  site  are  -33°  and  -40°  C, 
which  correspond  to  8.2  and  9.2  km,  respectively.  Both  altitudes  are  within  the  radar-derived  cloud  boundaries 
but  near  the  cloud  top  at  the  particular  WB-57F  overpass  time. 
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Figure  26  (a)  Partial  sounding  of  the  middle  and  upper  troposphere  obtained  using 
HARP  flight  above  Hanscom  AFB  (squares  correspond  to  retrieved  cloud 
temperature),  (b)  Comparison  of  cloud  top  heights  (obtained  from  radar/lidar 
measurements)  and  retrieved  cloud  temperatures.  Retrieved  temperatures 
are  displayed  for  Tracks  1, 3,  5,  and  8  (in  chronological  order). 


To  verify  the  retrieved  cirrus  cloud  optical  depths  and  ice  crystal  mean  effective  sizes,  we  need  to  use 
the  in-situ  ice  crystal  size  distributions  obtained  by  the  optical  probes  on  board  HARP.  Between  2030  and 
2040  UTC  (case  1),  and  between  2147  and  2157  UTC  (case  2),  HARP  flew  a  race  track  pattern  which  was 
within  30  km  from  Hanscom  Air  Force  Base  (Figure  27),  and  are  approximately  collocated  and  coincident  with 
the  Tracks  1  and  3  of  the  ARES  flights.  For  both  cases,  the  flight  altitude  was  at  about  8  km,  corresponding 
to  the  radar-detected  cirrus  cloud  altitude.  Ice  crystal  size  distributions  sampled  by  the  2D-C  probe  on  board 
HARP  for  both  cases  have  been  shown  in  Figure  6.  From  scattering  and  absorption  calculations,  visible 
extinction  coefficients  of  0.54  and  0.57  km'1  were  obtained  for  cases  1  and  2,  respectively.  Moreover,  based 
on  the  cloud  boundary  altitudes  derived  from  the  8.6-mm  radar  measurements,  the  averaged  cirrus  cloud 
thicknesses  during  the  two  periods  were  about  1.75  and  2  km  for  cases  1  and  2,  respectively.  Because  the 
vertical  profile  of  the  2D-C  median  volume  ice  crystal  diameter  obtained  during  the  ascending  run  of  HARP 
inHiratre  that  ice  crystal  sizes  vary  by  no  more  than  20  pm  between  6.5  and  8  km,  it  appears  reasonable  to  use 
the  computed  extinction  coefficients  based  on  single-level  aircraft  measurements  to  represent  the  vertically 
averaged  cirrus  cloud  extinction  coefficients.  The  computed  optical  depths  based  on  the  preceding 
approximation  are  0.94  and  1.14  for  cases  1  and  2,  respectively.  The  mean  retrieved  optical  depths  are  0.89 
and  1.42  for  Tracks  1  and  3  of  the  ARES  flight,  respectively.  The  retrieved  mean  optical  depth  for  Track  1 
agrees  relatively  well  with  that  from  the  2D-C  probe.  However,  for  Track  3,  the  difference  between  the  two 
optical  depths  is  on  the  order  of  0.3.  There  are  two  reasons  for  such  discrepancy.  First,  estimate  of  the  optical 
depth  from  the  2D-C  probe  is  based  on  the  assumed  composite  ice  crystal  size  distribution  at  a  single  flight 
level  within  the  cirrus  cloud.  Second,  the  variation  range  of  the  retrieved  optical  depth  for  Track  3  is  much 
larger  than  that  for  Track  1 . 

Finally,  we  use  the  computed  mean  effective  ice  crystal  size  based  on  single-level  aircraft 
measurements  to  represent  the  vertically  averaged  cirrus  cloud  mean  effective  size  and  utilize  a  cirrus  cloud 
model  composed  of  50%  bullet  rosettes,  30%  hollow  columns,  and  20%  plates.  For  each  shape,  a  different 
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relationship  between  the  maximum  dimension  (L)  and  width  (W)  is  employed:  for  plates,  W  2.02,  L  0.449, 
for  hollow  columns,  W  =  11.3,  L  0.414  for  L  >  200  pm,  W  =  8.479  +  1.002  L  0.00234  for  L  <  200  pm;  and 
for  bullet  rosettes,  W  =  2.3 103  L  0.63.  The  mean  effective  sizes  for  the  two  size  distributions  are  based  on 
the  average  of  the  mean  effective  size  for  each  shape  weighted  by  the  percentage  distribution.  They  are  57 .2 
pm  for  case  1  and  55.2  pm  for  case  2.  The  retrieved  mean  effective  sizes  are  64.2  pm  for  Track  1  and  53.3 
pm  for  Track  3.  Both  compare  reasonably  well  with  the  values  derived  from  the  2D-C  probe. 

4.  SUMMARY 

The  single-scattering  properties,  including  phase  function,  single-scattering  albedo,  and  extinction 
coefficient  for  representative  ice  crystal  size  distributions  in  the  spectral  range  2-6.4  pm  have  been  calculated 
using  the  unified  theory  for  light  scattering  by  ice  crystals  developed  by  our  research  group.  We  have 
investigated  the  range  of  ice  crystal  bounding  parameters,  including  the  refractive  indices  for  ice  in  the  ARES 
spectral  interval  and  the  ice  crystal  size  and  shape  for  midlatitude  cirrus  cloud  systems.  The  calculations  have 
been  carried  out  for  three  representative  cirrus  cloud  ice  crystal  size  distributions  corresponding  to  typical  cold, 
intermediate  and  warm  cirrus  clouds  in  the  midlatitudes.  We  have  also  investigated  the  effects  of  incorporating 
small  ice  crystals  on  the  phase  function,  single-scattering  albedo,  and  extinction  coefficient. 

Further,  we  have  developed  a  remote  sensing  algorithm  and  the  associated  computer  coding  (see 
Appendix)  for  the  retrieval  of  cirrus  cloud  temperature  (height),  optical  depth,  and  ice  crystal  size  based  on  the 
radiative  properties  of  cirrus  clouds.  The  ARES  5. 1-5.3  and  3.7  pm  channels  were  used  for  this  purpose.  This 
scheme  has  been  applied  to  the  ARES  data  collected  on  September  16,1995,  over  the  Hanscom  AFB  area.  The 
retrieved  cloud  temperatures  compare  reasonably  well  with  the  collocated  and  coincident  values  determined 
from  the  ground-based  8.6-mm  radar.  The  retrieved  optical  depth  and  mean  effective  size  also  match  the  values 
determined  from  in-situ  microphysical  measurements. 
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